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We consider here a manifestly covariant quark model of the nucleon and the A, where one quark 
is off-shell and the other two quarks form an on-shell diquark pair. Using this model, we have shown 
previously that the nucleon form factors and the dominant form factor for the 7A'^ — > A transition 
(the magnetic dipole (Ml) form factor) can be well described by nucleon and A wave functions with 
S-state components only. In this paper we show that non-vanishing results for the small electric (E2) 
and Coulomb (C2) quadrupole form factors can be obtained if D-state components are added to the 
A valence quark wave function. We present a covariant definition of these components and compute 
their contributions to the form factors. We find that these components cannot, by themselves, 
describe the data. Explicit pion cloud contributions must also be added and these contributions 
dominate both the E2 and the C2 form factors. By parametrizing the pion cloud contribution for 
the transition electric and Coulomb form factors in terms of the neutron electric form factor, we 
estimate that the contributions of the A D-state coupled to quark core spin of 3/2 is of the order 
of 1%, and the contributions of the A D-state coupled to quark core spin 1/2 is of the order of 4%. 



I. INTRODUCTION 

Understanding the internal structure of the baryons is 
both an experimental and a theoretical challenge. Ex- 
perimentally, the main source of information has been 
the electro- and photo- excitation of the nucleon, which 
allows to parametrize the baryons internal structure in 
terms of their electromagnetic form factors. Very accu- 
rate Jlab data [l|, [2, Q exist nowadays for the nucleon 
elastic form factors. Also, theoretical models for the nu- 
cleon form factors are able to describe this data well 
0, H, H, 0, S 111- The next step is the description of 
the nucleon excitations, starting with the A resonance. 
In the recent yea rs new precise data have been collected 
from MAMI [13, [H, LEGS [111, MIT-Bates [H and 
Jlab [11, [3 in the region < 6 GeV^ (g^ = _q2 jg ^-^^ 
squared momentum transfer). The A'^A electromagnetic 
transition (7A^ A) has a simple interpretation in terms 
of the valence quark structure: the A results from a spin 
flip of a single quark in the nucleon. It is then under- 
standable that the magnetic dipole multipole Ml domi- 
nates the transition for low , and that the electric E2 
and the Coulomb C2 quadrupoles give only small contri- 
butions, of the order of few percent. For large how- 
ever, according to perturbative QCD (pQCD) [3 flTf - 
equally important contributions from Ml and E2 are to 
be expected, but the scale for the outset of that regime 
is not yet known exactly. 

Several theoretical descriptions have been proposed for 
low, intermediate, as well as for the large transfer mo- 
mentum regions. These descriptions involve two in- 
gredients: the valence quark and the non-valence de- 
grees of freedom. The non-valence degrees of freedom 



are essentially the sea quark contributions which repre- 
sent quark-antiquark states, and are usually called me- 
son cloud effects. Due to its pseudoscalar character and 
its low mass, chiral symmetry assigns a special role to 
the pion [H, [l^, and pion cloud effects are therefore 
expected to contribute significantly to the baryon excita- 
tions. At low momentum transfer Effective Field Theory 
models based on chiral symmetry and perturbation the- 
ory (xPT) [13, HU, m, lil], with nucleon, A and pion 
degrees of freedom, and no internal structure consid- 
ered, work well. Effective Field Theories (EFT) describe 
the pion cloud effects at low momenta, but have a lim- 
ited range of apphcatiorK < 0.25 GeV^. At low Q^, 
the large limit 24, ,2^ can be used to establish the 
main dependence of the form factors and derive re- 
lations between the nucleon and the A'^A form factors 
[13. At large models within pQCD [H, ^ with 
quarks and gluons as degrees of the freedom, can be ap- 
plied. As for the intermediate momentum region, it may 
be appropriately featured by constituent quarks models 

M 



, D, 1281, |29|, 130, 131|, \M \M IM 1351, 1361, 133, IM 1391, |40|, 

and models based on hadronic degrees of freedorn, as the 
so-called dynamical models (DM) [U HI, jH, Hi, [ll, . 



Quark models with mixed coupling with pion fields have 


also been proposed [19, 27, 47, 48, 49, 50, 


51, 


52, 


53, 


5^. 



In the intermediate regime results from Vector Meson 
Dominance models j55{. QCD sum rules [s^ [stI and 
Global Parton Distributions (GDPs) [H, l59l. l60l|. have 
been presented as well. Finally, precise calculations 
are recently emerging from Lattice QCD calculations 
[bH . [6^ . m, [ill- For a review of the state-of-the-art in 
experiments and theory see Refs. [H, [gl, [66l. [67|. 

There is at present a strong motivation to pursue an in- 
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terplay between dynamical models and constituent quark 
models [l^ [13, O, [6^ . On one side, constituent quark 
models underestimate the result of the transition form 
factors, when not combined with explicit pion degrees 
of freedom [13, H IH, IH, Hi ■ On the other hand, dy- 
namical models are based on sets of equations coupling 
electromagnetic excitations to meson (de)excitations of 
the baryons, and include pion cloud effects naturally and 
non-perturbatively. Examples are the Sato and Lee (SL) 
[mils, [13, the Dubna-Mainz-Taipei (DMT) ^,'45^ and 
the (Dynamical) Utrecht-Ohio 46] models. Although 
very successful in the description of the the 7iV A 
form factors, they need to assume an initial parametriza- 
tion for the baryon transition vertex, interpreted as the 
bare vertex, where no pion loop is taken. In a less ad- 
hoc fashion, the bare vertex should therefore be derived 
from a quark model (66j . At the same time we can use 
the available models to extract the bare vertex. This is 
the goal of the EBAC program [H, [H, [6i|. Alterna- 
tive descriptions of the pion cloud and its relation to the 
deformations of the baryons were also proposed by Buch- 
mann et. al. [H, [tO, UM, [zl] • For an updated review of 
the dynamical models see Refs. [43[, [tsLIT^I . 

From the literature it is not clear which effects are due 
to the valence quarks, and which are related with the pion 
cloud, in particular for the E2 and C2 multipoles. There 
is a disagreement about those effects, between models 
based on different formalisms, as the DM and EFT mod- 
els, and even between models based on the same frame- 
work, as Effective Fields Theories [tI] . Also in the exper- 
imental sector there are some ambiguities. The form fac- 
tors are extracted using multipole analysis based on Uni 
tary Isobar Models (UIM) like MAID 
[Zl [Zi, [13, [Mi or Jlab/Yereven [H, 
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mUM, SAID 
each leading 

to different results due to the differences in parametriza- 
tions of the background and resonance structures, even 
considering dynamical models [ill, [H, [H, [3, [H, [21] in- 
stead of UIM [H, [H, [H [131. The ambiguities involved 
in the interpretations of the data are well illustrated by 
the differences between the CLAS results [l^ and the 
MAID analysis of the same data [t^, and also the re- 
cent preliminary CLAS data analysis [43, 67, 84] based 
on different models. The results that we obtained here 
for in particular, illustrate well the need to clarify 
these issues, as we will discuss later in this paper. 

In a previous work [1, [^, [6^ the spectator formalism 
[ssl . [s^ [STj was applied to the nucleon and to the A 
baryons, considering only S-state wave functions. As 
shown in that work, with S-waves alone in the baryon 
wave functions, only the dominant of the three form fac- 
tors for the — > A transition does not vanish. There- 
fore here we explore for the first time the effects of the D- 
states in the A wave function within that formalism, and 
show here that those components in the A wave function 
lead to non-zero contributions for the subleading form 
factors E2 and C2. The origin of the D-wave states is 
well-known: in the pioneering work of Isgur-Karl [28l [65| 
the baryons are described as a system of confined quarks. 



where a tensor color hyperfine interaction is generated 
by one-gluon-exchange processes. This tensor interac- 
tion leads to SU(6) symmetry breaking, and allows the 
transition from the ground S-state to an excited D-state. 

For each of the three NA electromagnetic transition 
form factors, we identified and separated the roles from 
the different partial wave components. While the mag- 
netic dipole form factor G^^, the dominant contribution, 
is mainly due to the transition between the nucleon and 
the S-state of the A, the electric quadrupole form factor 
G*^ proceeds through the transition to a D-state of the 
A corresponding to a three-quark core spin of 3/2. Fi- 
nally, the Coulomb quadrupole form factor becomes 
non-zero, only when the transition to a D-state of the A 
corresponding to a three quark core spin 1/2 is switched 
on. Nevertheless, and in agreement with other quark 
models, we conclude that the valence quark effects are 
not sufficient to describe the E2 and C2 data [H, [H [13] ■ 

Additional mechanisms involving the sea quark states, 
mainly the pion cloud effects, are needed to fill the gap 
between the theory and the experimental data. The sys- 
tematic and consistent treatment of the pion cloud mech- 
anisms is out of the scope of this work, which is focused 
on the D-state effects, but is planned for a future work. In 
order to estimate the magnitude of the D-states we con- 
sidered the simple parametrization of the pion cloud in 
terms of the nucleon (neutron) electric form factor, with 
no additional parameters. This parametrization was de- 
rived from the basic properties of the quark models (large 
Nc limit and also SU(6) symmetry breaking) and is lim- 
ited in its range of application to low Q^. Nevertheless, 
we need to include a pion cloud parametrization for a re- 
alistic estimate of the weight of the D-wave components 
in the A wave function. 

This paper is organized as follows: the formalism for 
the D-wave components of the A wave function is ex- 
plained in Sec. [TTl the definitions of the form factors and 
other general results are introduced in Sec. IIIIl the issue 
of gauge invariance and how it couples to the orthogo- 
nality of the initial and final state is dealt in Sec. IIV[ the 
formulas for the form factors within the valence quark 
model used here are given in Sec. [V] In Sec. IVII we dis- 
cuss the the contributions of the sea quarks (pion cloud), 
and in Sec. IVIII we present the numerical results for rep- 
resentative models based on valence and sea quarks. A 
discussion follows in Sec. IVIIIi and conclusions are pre- 
sented in Sec. [TXl 



II. NUCLEON AND A WAVE FUNCTIONS 

In the framework of the spectator theory [H, [s^ a 
baryon with four-momentum P is taken to be a bound 
state of a quark-diquark system, with relative four- 
momentum fc, and is described by a covariant amplitude 
^'(P, k). The diquark is taken to be on-mass-shell with an 
average mass rus- The 3-quark wave function amplitude 
depicted in the diagram of Fig. [T] has S-wave components 
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FIG. 1: Baryon quark-diquark wave function amplitude. 



which, for the nucleon and the A, were ah'eady presented 
in Refs. 0, [sl]. Therefore they wiU be only very briefly 
reviewed here, where the main focus is on the construc- 
tion of D-wave components within the same underlying 
formalism. In the following we will use H to denote ei- 
ther the nucleon (TV), with mass rriH — rriN = m, or the 
A with mass ttih = toa = M. 

The antisymmetry for the color part of the baryon 
wave function implies that, for S and D-waves, the spin- 
isospin part of the wave function is symmetric. This in 
turn implies that the diquark has positive parity. 



S-wave components of the Nucleon and A wave 
functions 



The S-wave part of the nucleon wave function has two 
components corresponding, respectively, to a diquark of 
spin 0-isospin and a diquark of spin 1-isospin 1. Label- 
ing the polarization of the spin-1 diquark by A, these two 
terms for the nucleon amplitude shown in Fig. [T] can be 
written 



N \„ 



(P,fc) 



1 



V2 



[(/.?u„(P,A„)-0}£^|,{/„(P,A„)] 



(2.1) 



The isospin states (f)j'' (with Id = 0,1 the diquark isospin) 
are respectively — and 0} = — -^(r • ^^*)x/ 

where ^'^* is the diquark isospin-0 state and ^i™* are the 
cartesian components of the isospin-1 state with projec- 
tions m = 0, ±1, and xi is the nucleon isospin state with 
nucleon isospin projection / = ±1/2. As explained in 
Refs. {t-^^*)xi generates the 3-quark isospin state 

in terms of the nucleon isospin. 

On the Dirac space, the the spin-0 component is sim- 
ply u„(P, A„) (denoted simply by u(P, A„) in our pre- 
vious work), where A„ is the projection of the nucleon 
spin along the z axis. The spin-1 component is a vector 
product of the diquark polarization vectors e"p and the 
Dirac operator 



u„(P,A„) (2.2) 



with (generalizing the definition of u„ to for an arbi- 
trary hadron with mass mn) 



2mH 



1 

2XhP 



where P„(P) 



(2.3) 

Note that Ujj{P,Xh) is 



independent of mn when P = 0. 

Note also that Ua{P,Xn) satisfies the two auxihary 
conditions 



(/'-m)C/,(P,A„) =0 
P"C/„(P,A„) =0. 



(2.4) 



The S-wave component of the A wave function must be 
symmetric in spin and isospin. Because the total isospin 
of the A is 3/2 the diquark spin 0-isospin component 
cannot contribute, and the A wave function can contain 
only a diquark spin-1 isospin-1 component. As defined in 
Ref. 68] the S-wave component of the A wave function 
is written 



vl/i(P, k) = -iPiiP, k)4>]ef*pWf,iP, Aa) 



(2.5) 



In this expression W/3(P, Aa) is the Rarita-Schwinger 
vector-spinor [ssl . [s^ satisfying the usual auxiliary con- 
ditions 



p0 



W/3 



0, 



I'^wp = , 



and the Dirac equation 

{p-M)wf,{P,XA) 



0. 



(2.6) 



(2.7) 



The function (pj — (T ■ ^^*)xi is the isospin part of the 
state, with T the isospin transition operator, and xi the 
isospin-3/2 state of projection /. 

For a particle of mass mn and three-momentum P in 
the z direction, the fixed-axis basis states used in (|2.ip 
and (12.51) are defined as 



j^P = T^(0,l,±*,0) 



r^P = (P,0,0, 

niH 



''H 



p2). 



(2.8) 



These are the diquark fixed-axis polarization states dis- 
cussed in great detail in Here it is sufficient to note 
that they satisfy the orthogonality condition P -exp = 0, 
and that 



pii pi' 



(2.9) 



Because of the orthogonality condition, the wave func- 
tions (|2.ip and (|2.5p satisfy the Dirac equation for the 
mass niH- 
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As discussed in "s"] the fixed-axis diquark polarization 
states introduce no angular dependence in the wave func- 
tion, and therefore are convenient to describe S-states, 
without introducing any extra constraint. We will show 
subsequently here that they are also convenient to define 
higher angular momentum states. 

Finally, for the scalar wave functions ipff and tp^ in 
(|2.ip and (|2.5p , which describe the relative quark-diquark 
radial motion, we use the parametrizations 



where 



XH 



No 


ms{(3i + 


■ Xn){P2 - 


"tXn) 




Ns 




ms(ai 4 


XA){a2 


+ Xa)2 


[P-kf 


- [ttih - 


m^)2 



mnms 



(2.10) 
(2.11) 

(2.12) 



for H = N, A. This parametrization allows for an in- 
terplay between two different momentum scales in the 
problem, quantified by the /3i and P2 parameters for the 
nucleon, and ai and a2 for the A. 

The factors A^o a-nd Ns are normalization constants 
fixed by the condition 



/ [V^(P,fc)]' = l, 

Jk 



(2.13) 



with P = (m//, 0, 0, 0) and the covariant integral is 
defined in Eq. (|4.3p below. The normalization condition 
is consistent with the charge constraint 

Qz = 3^ / *|(P,fc)ji(0)7°*i(P,fc), (2.14) 

where ji{0) is the isospin part of the quark charge op- 
erator, and Qi is the charge of the state with isospin 
projection / for either the nucleon {H = N with isospin 
1/2) or the A {H = A with isospin 3/2). See Refs. 
for more details. 

In this paper the A, with total angular momentum 
3/2, is composed of two positive parity subsystems (the 
spin 1/2 quark and the spin 1 diquark). We will refer 
to the total spin S of the state as the "core" spin, in 
order to distinguish if from the total angular momentum 
of the state (also called the "spin" of the particle). If 
the orbital angular momentum between the quark and 
diquark is zero, then the "core" spin of the A must be 
3/2. However, a state of positive parity and total angular 
momentum 3/2 can also be constructed if the orbital an- 
gular momentum of the constituents is a D-wave (L = 2) 
and the core spin is either S = 1/2 or 5 = 3/2. (For the 
nucleon, in contrast, the L — 2 orbital state can couple 
only to the core spin S — 3/2.) The next subsections will 
define these two possible A D-states. 



B. The two different spin-core D-wave components 
of the A wave function 

1. D-wave operator 

Turning to the construction of the momentum-space 
part of the D-wave components of the wave function of 
the A, we start by noting that, in relativity, inner prod- 
ucts of three vectors (and, consequently, magnitudes of 
angles) are not Poincare invariant. Hence, the operator 
T)"^ that generates a D-wave in the relative momentum 
variable k, only has the pure D-wave structure in the 
baryon rest frame. In any other (moving) frame the in- 
trinsic D-state will generate components in other partial 
waves. Therefore, we start by defining that operator in 
the rest frame of the A. To find its form in that frame we 
exploit the two features that define a D-wave: (i) V"'^ is 
bilinear in the 3-momentum vector part k in the A rest 
frame, and (ii) the integral of V"^ over all the possible 
directions of k has to vanish in the rest frame. 

It is convenient to introduce a four- vector that reduces 
to the 3-momentum k in the A rest frame. Defining this 
vector for an arbitrary hadron 



P ■ k 

k"" = k'^ —P" 



(2.15) 



where, in the hadron rest frame, k = (0, k) and fc^ — — k^. 
In terms of this vector, the two defining properties of V 
lead immediately to the operator 



V"f^{P,k) = ~k"k'^ - —g"^, 



where 



pa p/3 



(2.16) 



(2.17) 



(In this discussion we suppress the subscript H on both V 
and fc, relying on the reader to infer the correct operator 
from the context.) Note the constraint conditions 



(2.18) 



It is convenient to work with the spherical components 
of V, defined to be 



(2.19) 



Using the definition (|2.8p of the fixed-axis polarization 
states, it is easy to see that D is a hermitian matrix. 
In the hadron rest frame the matrix elements of D are 
related directly to the spherical harmonics l^j(k) with 
L = 2. A representation convenient for later applications 
is 



k2y4(k)(lA2m,|lA') 



(2.20) 
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where the vector couplmg (or Clebsch-Gordon - GC) co- 
efScients are denoted 

O'lMi i2M2|ji2Mi +M2) = C{jij2ji2; ^J.l^J.2) ■ (2.21) 



Equation (|2.20|) shows how the operator D can be inter- 
preted as the projection of the incoming direct product 
state of orbital angular momentum L = 2 (g) a spin-1 
vector , onto an outgoing vector state e*y 



AA' 



^ el, (1 A 2 m, I 1 A' ) if,^ (k) . (2.22) 



AA' 



2. Spin Projection operators 

To prepare for that construction of the D-wave compo- 
nents of the wave function, we recall the definitions of the 
spin-projection operators V1/2 and P3/2 previously used 
in Ref. [6^ (and in other works). These are constructed 
from the operator (|2.17p and the operator 



T) pa 

1 ==7 -— 2- 



•■H 



The operator has the property 



(2.23) 



(2.24) 



In terms of these operators, the projection operators can 
be written 



1 



(^1/2) 

(^3/2)"^ = - {Vi/2r^. 



(2.25) 



For details see Refs. [gSj, |90|, l9l| . We note that 7^3/2 can 
be cast into the form usually found in the literature. 



V(/'7"-P'^ + ^"7^ Z*)' (2.26) 



3771 



and that these spin projectors satisfy the closure and 
orthogonality relations 



(^1/2)"^ + (^3/2)"'' = 



{V^/2r^{V^/2)p^ - - (7'3/2)"''(^l/2)^^- (2.27) 



af}. 



Denoting the operators (|2.25p generically by Vs^ one eas- 
ily sees that 



PaiVsT" = 

7« ^3/2 = = ^3/2 7/3. 



(2.28) 
(2.29) 
(2.30) 



Note also that the state functions previously introduced 
in Eqs. (|2.1[) and (|2.5|) satisfy the expected eigenvector 
equations 



(^l/2)/f// 
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iV3/2)Jwp 



(2.31) 



Construction of the two possible D-state components of 
the A 



Using the operator V and the spin projection operators 
Vs introduced above, we can now construct D-state wave 
functions for the A. Just as the S-state wave function is 
a matrix element of the A initial state with a final state 
consisting of a quark and a diquark in a relative S-state, 
the D-state wave functions are matrix elements of the A 
initial state with a final state consisting of a quark and 
diquark in a relative D-state. The construction is carried 
out in two steps. First a D-wave dependence is intro- 
duced by contracting the V operator with the elemen- 
tary S-wave Rarita-Schwinger wave function Wa, giving 
the state 



(2.32) 



The resulting state W" satisfies the Dirac equation. 
Next, using PaV"^ = and the completeness relation 
Eq. (|2.27p we conclude that this is actually the sum 
(only) of two independent spin components 



0' 
{V1/2) 



';3 + (^3/2)"/3 Wf^(F,fc).(2.33) 



This leads to the definition of two independent D-wave 
A wave functions 

cl^D 25(AAa) = ~3e{; {Vs)pa WxJP^k) 

= -3e^; {Vs)0a V^^{P, k) w^{P, Aa) , 

(2.34) 

where the factor of —3 has been added for convenience, 
5=1/2 or 3/2, and e\ describes the state of the outgoing 
diquark, just as in the S-state wave functions Eqs. (|2.ip 
and (|2.5p . Equation (|2.34p defines the spin part of the 
two D-state wave functions only; isospin and radial parts 
will be added below. These wave functions satisfy the 
Dirac equation (|2.7p . 

It is interesting to see how the wave functions (|2.34p 
have the correct spin structure corresponding to the two 
different {L,S) couphng configurations, (2, i) and (2, |), 
both giving total J = 3/2. Here we summarize the main 
points; details are given in Appendix [Xl The first step is 
to introduce core spin wave functions (direct products of 
the spin-1 diquark and a spin-1/2 quark) with 5 = 1/2 
or 3/2. These core wave functions, denoted generically 
by Vs, are constructed using CO coefficients 



Vs"iP,Xs 



A 



(iAlA'|5A,)e^,pMA(P,A), (2.35) 
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where ua was defined in Eq. (|2.3p (with H ^ A). It 
is easy to see that the Vg satisfy the Dirac equation 
(|2.7p . that PaVg = 0, and it can be shown that they 
are eigenstates of the projectors Vs- The fact that V3/2 
also satisfies the special spin 3/2 constraint •jaVy-) — 
is shown in Appendix [X] 

These wave functions are orthonormal and complete 



Vs{P,X)Vsa{P,X') = Sxy 



As 



mA+ P 
2mA 



(2.36) 



Using the Dirac equation to introduce the projection op- 
erator into Eq. (|2.34p . and then inserting these expan- 
sions, allows us to express the D-state wave functions in 
the form 

<t>D2s{X) = -ie1*pY.^s piP,Xs) 

As 



{Vsa{P,Xs)'D''^W^{P,XA)} 



= (-1) 



s-i 



(2 mi-SXs\lXA)Yl^ysp{P.Xs)- 

(2.37) 



This displays the two states as sums over either an 
S —1/2 or 3/2 core wave function Vs times an orbital 
angular momentum L = 2 spherical harmonic function 
coupled to a spin 3/2 A state, and is demonstrated 
in Appendix |^ Using (|2.37p and the normalization of 
the Vs states, the normalization of the 4>D2S are 



1 

47r 



A 



|0i?2s(AAa)|' = k' 



^^(2 toz;5A,||Aa) 



(2.38) 



Using the definitions (|2.34p and adding the isospin fac- 
tor and scalar wave function, the complete A D-state 
wave functions are 



^'i\P,k) ^ ^ni{X) 4>] ^KHP,k) 
^^^{P,k) = 4,D3{X) 4>] KHP,k), 



(2.39) 



where the following simple forms were used for the scalar 
functions ip-^^ and V'a'^ 



V'fi(P,fc) = Nni 



1 



"^s(a3 + xa) 
Xdi 



Nd3 



Xa) 



(2.40) 
(2.41) 



The Dl state has two range parameters (as and a^) and 
the D3 state only one (as). The three parameters are 
adjusted to the data. We anticipate that our numerical 



results for the range parameters are consistent with an 
expected longer range (in r-space) for the D-states rela- 
tively to the S-state. Note that the definition of V guar- 
antees also that the D-state wave function will go as 
when k — > 0, as expected [o^l, and an additional mass fac- 
tor mj'^ is introduced to compensate for the dimensions 
introduced by this k^ dependence of the V matrix (so 
that the product of T> with or ip^^ have no dimen- 
sions). The power 4 in the denominators of the previous 
equations was chosen to reproduce the expected pQCD 
behavior for large (^^ _ i/q4^ q*^ _ ^/q6^ 

and also to assure the convergence of the normalization 
integrals 

Combining Eqs. (|23)) and (|2?39l) . the total A wave 
function can be written 



(2.42) 



where a and b are admixture coefficients. The Dl com- 
ponent with core spin 1/2 is orthogonal to both of the 
other components because of the orthogonality condition 
(|2.27p and the two components with core spin 3/2, ^f^"^ 
and ^f^, are orthogonal because the overlap integral is 
linear in /j.l2o(z) — 0. We will chose to normalize the 
individual states to unity, giving N = l/Vl + a2 + 62 for 
the overall normalization factor. 



C. Normalization and orthogonality condition 

The individual D-wave scalar wave function will be 
chosen to satisfy the normalization conditions 



= 1. 



(2.43) 



This determines the coefficients A'^^ii (as a function of 
Xdi) and Nds- 

The two components and are orthogonal to 
the nucleon S-state '^ff, but the component is, in 
general, not orthogonal to the nucleon S-state. This hap- 
pens because both wave functions have a core spin S = 
1/2, and even though the Dl state depends on Yj^^, it is 
impossible for both particles to sit simultaneously in their 
rest frame, so the angular integral always has some other 
angular dependence that prevents it from being exactly 
zero. The orthogonality condition 

J2 I ^^HP+,k)^NiP^,k) = 0, (2.44) 



where P+, P_ represent the baryon momenta for — 0, 
must be imposed numerically, and this can be done only 
at one value of Q^. This condition determines Xjji . As we 
will see below, our treatment of gauge invariance requires 
that we impose the condition (|2.44p at the point = 0. 
Working in the A rest frame, the momenta P+ and P_ 
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are therefore 



P+ = (M, 0,0,0) 



M2 



m 



2M 



0,0, 



M2 



2M 



where P+ (P_ ) is the momentum of the A (nuclcon) , /' 
(/) the isospin projection of the A (nucleon), and the 
operator T^'^ can be written in general jl^l as 

(2.45) r'3^(p,g) =GigV + G2/P^ + G'3/g^-G4/^. (3.2) 



To determine the coefficients Ndi and A/ji in the Df 
component, we first fix and a^. Then A^ii is deter- 
mined by (|2.44|) . and finally the value of Ndi fixed by 
the normalization condition (|2.43p . 



D. 



Properties of the wave functions under a 
Lorentz transformation 



The form for the wave functions given in Eqs. (|2.39p 
holds only for the case where the particle is mov- 
ing along the z direction (with 4-momentum P ~ 

(y^m|^~-t-I^, 0, 0, P). The generic wave function can be 
obtained from an arbitrary Lorentz transformation A: 

P'^ = Af'^P". (2.46) 

Under a Lorentz transformation we obtain: 

w'p{P') = k0"S{K)w^{P) 
u'{P') = S{K)u{P) 
V^p{P\k')^KKPV„p{P,k) 

S'\A){V'sUS{A) = A-A;(P5).p, (2.47) 

where u' and represents the states in the arbitrary 
frame. For simplicity, the dependence of the spinor states 
on the Wigner rotations acting on the polarization vec- 
tors has not been shown explicitly, and (Ps) are the 
projectors of (|2.25p with {P'g) the same projectors with 
P' = AP, one obtains the transformation law 



Z'JP',k')^S{A)A0^Z^{P,k) 



(2.48) 



for any vector spinor state Z. Finally, from (|2.48p the 
transformation laws for the total A wave function follows 



^'^{P\k')^S{A)^A{P,k). 



(2.49) 



In conclusion, we may derive the baryon wave function 
in any frame, where the four-momentum P is arbitrary, 
by means of a Lorentz transformation A on the wave 
function defined in the baryon rest frame. 



III. FORM FACTORS FOR THE -/N 
TRANSITION 



A. Definitions 

The electromagnetic iVA transition current is 

J'' - Wp{P+)rf''^{P,qh5u{P-)S,„, (3.1) 



Although we have omitted the helicity indices for these 
states, the transition current depends on both the he- 
licities of the final and initial baryons and on the pho- 
ton helicity. The variables P and q are respectively the 
average of baryon momenta and the absorbed (photon) 
momentum: 



q = P+-P- 



P_ 



(3.3) 



The form factors Gi, i — 1 , . . , 4 are functions of 
= —q^ exclusively. Because of current conservation, 
Q/iT'^^ = 0, only three of the four form factors are inde- 
pendent. In particular, we can write G4 in terms of the 
other three form factors as 



G4 = [M + Tn)Gi + 



A'P - 



-G2 



(3.4) 



and adopt the structure originally proposed by Jones and 
Scadron [9^. Alternatively (see below) we can write G3 
in terms of the other three 



(M + m)Gi + 



-G2 — Ga 



(3.5) 



The paramctrization (j3.2p in terms of the form factors 
Gi is not the most convenient one for comparison with 
the experimental data. More convenient are the mag- 
netic dipole (M), electric quadrupole (E) and Coulomb 
quadrupole (C) form factors. These can be defined di- 
rectly in terms of helicity amplitudes [l^, [9^ . Note that 
the form factor G3 does not enter directly into the ex- 
pressions for the helicity amplitudes because e^*9p = 



for all A. But, if we use the constraint (|3.4p to eliminate 
G4, G3 appears in these expressions and we obtain 



[{3M + m){M + m) + 



Gi 

M 



[M^ - m^)G2 
(M^ - m2)G2 



2g'G3 } 



2\Gl 



^ ' M 



2Q2G3 } 
K |4A-f Gi + (3M2 +m^ + 
2(M2 - m2 - Q2)G3 I , 



Q')G2 



where 



3(M + m) ■ 



(3.6) 



(3.7) 



(3.8) 



(3.9) 
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These three form factors G* (a — M, E, C) are, respec- 
tively, the magnetic, electric and Coulomb (or scalar) 
multipole transition form factors. 

As G*jy,j dominates at low momentum Q^, the following 
ratios are useful 



conservation condition (|3.5p . the numerator (at = 0) 
is proportional to 



™2 

(M + m)Gi + G2 - G4 



Q^G3, (3.17) 



Rem{Q^) — — 



and 



i?SA/(Q' 



|q| G*c{Q^) 



(3.10) 



(3.11) 



2MGl,{Q^y 

where q is the photon 3-momentum in the A rest frame 



which, if G3 is finite, approaches zero as 0. 

We prefer the independent choice Gi, G2, G4 because 
it enables us to discuss the restrictions imposed by cur- 
rent conservation in a more transparent way. Many mod- 
els do not automatically conserve current (this is true for 
our Dl component, as we will discuss below). If we start 
with a model that does not naturally conserve current, 
we prefer to impose current conservation by modifying 
the current in the following way: 



with 



2M ' 



d± = [M ± mf + 



(3.12) 



(3.13) 



^ 



(3.18) 



The analysis of the transition at large in the pQCD 
regime (where quarks and gluons are the appropriate 



de grees of freedom) gives G^j 



and 



B. 



The GiyG-ZjGs set versus the Gi,G2,G4 set 



The representation of the electromagnetic current in 
terms of the 3 independent (Gi,G2,G3) form factors, as 
proposed by Jones and Scadron [93| , is not the most con- 
venient choice that can be made. As mentioned above, 
the form factor G3 is not part of the helicity transition 



This way of imposing current conservation is, of course, 
not unique, but has the nice property that the additional 
term added is proportional to q^, and hence makes no 
additional contributions to any observables obtained by 
contracting the current with another conserved current 
or with a photon polarization vector, always orthogonal 
to (in the Lorentz gauge, our choice). In the applica- 
tion discussed in this paper, the modification p.lSp will 
only alter the G3 form factor, and when we use the ex- 
pressions (|3.14p - (|3.16p we see that they are unchanged 
by any modification of G3. Hence, our method allows 
us to choose G3 to satisfy current conservation, without 
changing the basic predictions of the theory. 

However, current conservation is like the Cheshire cat, 
while the consequences of imposing it seem to have van- 
ished, a "smile" still remains. What remains is the re- 



amplitudes given by the operator e^(g)J'', due to the quirement that there is no singularity in G^ as ^ q. 



condition e ■ q = 0. For this reason it seems natural to re- 
place the set (Gi,G2,G3) by (Gi,G2,G4). In this basis. 



G*M: G^, and G^ are are given by 



M — 
G*E = 

Gh = 



n[2G, + d+^] 

.{2G4-d,^} 
^{2(M-m) d+Gi 
-2(1/2 - TO^ - ( 



(3.14) 
(3.15) 

- rf+rf_G2 

^)G4} (3.16) 



This requirement is satisfied by modifying the form fac- 
tors in such a way that the linear combination (|3.17p is 
zero at = 0. Implementation of this requirement will 
be discussed below. 



C. Simple relation for Gc 

In the following discussion we will work in the rest 
frame of the outgoing A, where the four- momenta 
become 



where d± were defined in Eq. (|3.13p . Note that the mul- 
tipole form factors G^j and G*^ do not depend on G2. 

Eq. p.l6p for G^ presents an apparent singularity 
when = 0. The presence of this apparent singular- 
ity is the historical reason for choosing Gi , G2 , G3 to be 
the independent form factors; this choice gives finite form 
factors under any circumstances. However, if the theory 
conserves current, with a G3 that is finite at = (a 
required feature of any consistent model), then the sin- 
gularity disappears as —> 0, since, using the current 



<z'' = (^,0,0,|q|) 

„ f2M-Lj Iql 

= 0- — 

I 2 ' ' ' 2 



(3.19) 



where |q| was given in Eq. (|3.12p . and the photon energy 
Lj can be written in terms of the nucleon energy oj — 
M - + |q|2, or 



P+ ■ q NP - m 



2 _2 /o2 



M 



2M 



(3.20) 
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In this frame the photon moves in the +z direction, with 
polarization vectors 



1 



=T-=(0,l,±i,0) 



1 

Q 



eo, = 7^(|q|,o,o,^). 



(3.21) 



Note that the transverse states (A = ±1) are identical 
to those defined in Eq. (|2.8|1 . but that the longitudinal 
state is very different. All of these satisfy the constraint 
g^e^ — 0, and because > are identical to helicity 
states. While we will work out the explicit relations in 
this rest frame, all relations that are derived from four- 
vector scalar products are, of course, independent of the 
frame. 

Introduce the photon helicity amplitudes of the electro- 
magnetic transition current (|3.ip (for a general discussion 
of helicity amplitudes see Refs. 



Jn — J\ 



(3.22) 



where the polarizations of the N and A will remain un- 
specified. Note immediately that Ga does not contribute 
to any of these amplitudes, and because P ■ e± =0, the 
transverse amplitudes do not depend on G2- The only 
amplitude that depends on G2 is the longitudinal JTq. 
Using the relations 

|q|M _ 1 
Q - ap 

e^,=a,,--HapP- (3-23) 

we can reduce the terms and e^^ that occur when 
using (13. ip to evaluate Jo, and obtain 



_ (M + m) 3Q 



(3.24) 



where s and s' are, respectively, the nucleon and A spin 
projections along the z-axis, and 



TZs's = u)0(P+,s')/75w(P-,s) 



Sss'{2s) 



2d+ d_ 
3mM4M' 



(3.25) 



We emphasize that, provided we use Eq. (|3.16p to define 
G^, this relation holds for all models, even those that do 
not conserve current. Note that TZ ^ only if the spin 
projections are equal (s — s'). We may conclude that 
Gq only if (for example) JTq ii 7^ 0, and using Eq. 
PTM)) we obtain 



Jo. 



(M- 



4m 



3 d- 
2mM 



M 



(3.26) 



IV. THE ELECTROMAGNETIC CURRENT 
WITHIN THE SPECTATOR MODEL 

In this section we study how the electromagnetic cur- 
rent can be constructed within the constituent quark 
model (CQM) for the baryon structure presented in 
Secini 

In any CQM model the quarks making up the baryons 
are not point particles, but composite valence quarks, 
dressed by their gluon and sea quark structure. Here 
we use the Covariant Spectator Theory and assume the 
baryon is a quark-diquark system, as explained in Sec. [Til 

The on-shell diquark mass rUg scales out from the elas- 
tic form factor, which turns out to be independent of the 
diquark mass [J| . This mass does not scale out of the deep 
inelastic results (DIS) and the qualitative description of 
DIS leads to the estimate of ~ 0.8to, allowing a nat- 
ural interplay between low and high energy phenomenol- 
ogy. This interplay is needed since the factorization into 
low and high energy scales does not apply exactly. 

In the following we will explain how gauge invariance 
conveniently constrains the current, when the internal 
structure of the quarks is parametrized in terms of phe- 
nomenologically fixed wave functions. 



A. Implications of the choice of current 

1. Simple current 

Constituent quarks are dressed particles with a com- 
plex effective structure, an effective charge and magnetic 
moment. Therefore, their current consists of a Dirac and 
a Pauli term, and can be written as 



J2- 



2m 



(4.1) 



(The subscript "a" on the current will be dropped in 
subsequent discussion, and will be used only when we 
need to distinguish this current from the modified current 
discussed in the next subsection.) The form factors ji 
and j2 are normalized in order to describe the nucleon 
charge and magnetic moments (as functions of the quark 
isospin /) as discussed in Ref. [4|. The explicit formulas 
are defined by the Eqs. (I4.15P and (|7.ip below. The quark 
current (|4.ip is not of the most general form. In the 
next subsection we will consider a more generic case, in 
light of the discussion on gauge invariance that unfolds 
immediately here as consequence of (|4.ip . 

To start this discussion, given the quark current (|4.ip 
and the nucleon {'^n) and A (^'a) wave functions, we 
write the transition current between these states. With 
a positive parity axial diquark the only allowed states 
for the nucleon and A are S and D states, since P-states 
are ruled out (unless they are associated with the lower 
relativistic components, not discussed so far in this series 
of papers). 
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the current (|4.2p as 



P. 



+ 



TV, A 



/ 



k 



TV 



FIG. 2: Relativistic impulse approximation. 



To simplify the formulas we will exclude the isospin 
from the discussion (later in this paper we show how to 
include the isospin explicitly) . In impulse approximation 
0, m, S S (H, the transition current takes the 
form 



■9 



pp. 



= 3j.V / ^Ar^N~ij2—p{Q'), (4.6) 



where 



M + m 



(4.7) 



These equations hold in any frame. 

Next, using Eq. (14. 6p . the relations (|3.23p . and the 
Dirac equation, we find an alternative form for the lon- 
gitudinal current 



\ J k 



(4.8) 



(4.2) 

where 



where all momenta and spin projections (s' for the A 
and s for the nucleon) have been suppressed. The factor 
3 sums up the contributions of the three quarks, the sum 
is over all intermediate polarizations A of the diquark, 
and 



(4.3) 



is the covariant integral with Eg = ml + as the di- 
quark on-mass-shell energy. The initial and final momen- 
tum dependence are not explicitly included for simplicity. 

As discussed in ^ equation (I4.2p . for the states and 
^'at defined here, goes beyond the scope of the relativistic 
impulse approximation (RIA) shown diagrammatically in 
Fig. [21 and includes some effective two body currents. 

Since both the final and initial states satisfy the Dirac 
equation: f+^A = M^a and /'_*Ar = m^N, the Pauh 
current can be simplified using the Gordon decomposi- 
tion 



' 2m 



M + m 
2m 



J2 / *A7''*w(4.4) 



{P+ + P- 
2m 



J2 I '^A* 



The last term is proportional to (inserting the spin pro- 
jections for clarity) 



Ps's{Q^) *As' '^Ns- 

\ J k 



(4.5) 



With this definition (dropping references to s' and s 
again), we can use the Gordon decomposition to write 



2m{M + m) 



(4.9) 



Both the current j7o and the factor p depend on the spin 
projections s' and s, suppressed so far. Taking the spin 
projections s = s' = ^ and combining this result with 
Eq. (|3.26p gives in terms of ^ 



AmM / &Mm 



(4.10) 



This result holds in any frame. 

Now we connect some of these results to the diver- 
gence of the simple current (|4.ip . Noting that the Pauli 
term, proportional to is automatically conserved, the 
divergence of current depends only on the behavior of 
the Dirac term, proportional to ji . Evaluating the diver- 
gence gives 



q ■ 



= 3{AI-m)n p{Q^)- 



(4.11) 



Because the masses are different {M 7^ m) this (frame 
independent) result shows that the simple current (|4.ip 
will be conserved for an electromagnetic transition from 
a state vPat to a state ^I^a if and only if p[Q^) ~ 0. We 
showed in Ref. [6^ that this term will vanish identically 
(for all values of Q^) if the core spins of the two states 
are different. This is true for the nucleon S-state to A 
S-state transition, and also for the transition from the 
nucleon S-state to the A D3 state. It is not true for the 
transition to the Dl state, as discussed briefly above, and 
in more detail below. 
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Combining Eqs. (|4.1ip and (|4.10p gives the following 
interesting connection 



2 2 



3(M 



Jc 



6mM AraM 



G. 



c- 



(4.12) 



The consequence of this equation is that the simple 
current (|4.ip will be conserved if and only if G*^ — 0. 
Alternatively, since transitions to the A S and D3 states 
conserve the current (j4.ip . these cannot give a non-zero 
G*(j. To build a model in which 7^ we must find a 
different current, and this leads us to the next subsection. 



definition of the current through the following isoscalar 
and isovector decomposition, as in Ref. [4]: 



(4.15) 



where i = 1,2 and /i± and /2± were adjusted by the 
charge and magnetic form factors of the nucleon and were 
normalized to /i±(0) = 1, /2±(0) = k±. Only the isovec- 
tor form factors, /i_ contribute to the ~* A transi- 
tions. 

The overall isospin factor can be calculated separately, 
and was worked out in Ref. 1681. This factor is 



2. Modified current 



Following a previous work [ 
rent ((i7T|) by 



J I = ]1 



r 



32- 



we replace the quark cur- 



j^, + Ajf. (4.13) 



2m 



It is easy to evaluate the additional term. 



^/3 , 



(4.16) 



All formulas derived in the precious sections still valid if 
we replace 



3i 



Go fi- 



(4.17) 



Ajj' = 3(M 



(4.14) 



This shows that all the good properties of the previous 
current remain intact; when p = for all values, 
the Dirac current Ji7'^ is conserved and the correction 
term vanishes identically. The advantage of the current 
(|4.13p is that current conservation is guaranteed, even 
when p{Q^) does not vanish identically. 

The only possible problem with the new current is that 
it might be singular at — 0. This singularity must be 
removed by imposing the requirement that p{Q^) — >■ 
as Q"^ 0. Since eo? • <? = 0, Eqs. ([48]) and ((iTTOll are 
unchanged, and this requirement also means that is 
finite at = 0, guaranteeing that the apparent singu- 
larity in Eq. p.l6p does indeed cancel. 

The condition that guarantees that p — at = 
was already introduced above in Sec. Ill C[ Eq. (|2.44p . 
The importance of this orthogonality condition was em- 
phasized in Ref. (98,] ; imposing it ensures that the current 
(I4.13P is well defined and conserved for all Q^. 

To summarize: in the present model the orbital angu- 
lar momentum states are not derived from an underlying 
hamiltonian. Therefore the A state with {L,S) — (2, i), 
with the same core spin quantum numbers as the nucleon 
state, even though carrying the correct spin-isospin sym- 
metries, does not have a spacial scalar part V'a^ ^^^^ 
ab-initio orthogonal to a nucleon state. The orthogonal- 
ity is imposed by a judicious choice of the parameter Aui 
in Eq. 1^:10)) . 



B. Isospin dependence of the current 

For simplicity, we did not include isospin in the dis- 
cussion in the previous subsection. It is included in the 



V. VALENCE QUARK CONTRIBUTION FOR 
THE FORM FACTORS 



The impulse approximation for jN A transitions 
from the nucleon S-state to each of the A states can be 
written, using Eq. (|2.42|) 



a J^3 



(5.1) 



where the index identifies the A state. From this we 
can calculate the form factors Gi, G2, and G4 defined in 
Eqs. (13. ip and (|3.2p . and using the definitions given in 
Eqs. (I3.14p - (|3.16p calculate the multipole transition form 



factors G\,j , G^ and G^ . 



A. Transition to the A S-state 

The transition current from the nucleon S-state for the 
A S-state was already evaluated in the Ref. [1^. Using 
the upper index S to indicate the A state, the results are 



3\/3 M + m 

gUQ') - 

Gg(Q2) ^ 0, 
where j„ is the analogue of (|4.7p 

_ M + m 

Jv — /l- H 7^ /2- 

Zm 

and 



fvis, (5.2) 
(5.3) 

(5.4) 



/ V^i(P+,fc)V^(P-,fc), (5.5) 
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is the overlap integral of the radial (scalar) wave func- 
tions. Asymptotically we have Gf^ ^ as showed 
in Ref. |68| . 

According to Eqs. ((XT^ and (PTTB)l . Gf and G§ van- 
ish because the terms involving Gi, G2, and G4. cancel 
exactly. 



This general expression may be reduced further by noting 
that, in a coUinear frame in which none of the momenta 
have components in the x or y directions, the only de- 
pendence of the integrand on the azimuthal angle (f is in 
the angular dependent term T). Hence, we may average 
over this angle using the (covariant) identity 



B. Transitions to the A D-states 



1 

2^ 



— / d^V''PiP+,k) = b{k,q)R 



.a/3 



(5.12) 



The transition currents to the D states are 

Jd2S = 3E / ^^''{P+,m^NiP^,k) 
\ •J h 



A 

= WMP+)r^'^25(F,g)75"„(P-)'5.... (5.6) 

where we suppress all reference to the spins of the nucleon 
and A. Substituting for ^fjv using Eq. (|2.ip . \E'£^'^ using 
Eqs. (|2.34p and (|2.39p . and using the general reduction 
(|4.6p gives 



(^-^"'^ ^Vf^"^^, (5.7) 



where A^" is the sum over the fixed axis diquark polar- 
izations (previously derived in Ref. [sl. Issj) 

A a' a \ ^ ^a' 



X 



pa pa 



pa' 



AP 



pa 







(5.8) 



with 



Mm 



a 
b 



b{Mm + b) 
P+P-. 



(5.9) 



and the two current operators emerging from the reduc- 
tion (14.611 are 



-h-- 



pfj, 



(5.10) 



Using the conditions (|T4)) . ((2^ . and ((2?28| . the part of 
the expression for T£)2S in curly brackets { } reduces to 



where 



KM) 4^ 4? 



2 



1 



~af3 



(5.13) 



with k and q defined as in Eq. (|2.15p [with the substitu- 
tions P P+ and mjj — > M] . This identity is proved in 
Appendix |B] 

Using the conditions ()2.6|) and (12.281) again, the (p av- 
erage of (|5.1ip can be simplified 



—^{Ps)P'a 



mM + b 



-5 ql^ 

2 "2s,3 y a 



(5.14) 

This will now be evaluated for the two cases of interest. 
1. Nucleon(S) A(D3) 



The term in round brackets in Eq. (|5.1ip commutes 
with O2 (an identity operator on the Dirac space), and 
hence, for the transition to the spin 3/2 core state (D3) 
with Vs — P3/2 gives zero (this is the p term discussed 
above). Commuting the term in round brackets through 
Oi, letting — > M when it operates to the left, gives 



7^ I T 



P^[P+ -M] 
mM + b 



p^[2P_i 



2M-^^■ 



mM + b 



(5.15) 



For the S=3/2 case under consideration, the 7^7^ terms 
vanishes, and combining this with the remaining terms 
gives 



1 



— (,P3/2)/3'Q - off a 

3 



X 7 



„ P^[p+-M] 



mM + b 



(5.11) 



X g 



P^[P^ - M7'' 

mM + b 



(5.16) 
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Now, we know that the terms proportional to 5^ can be 
ignored (they determine G3 which we already know is 
just right to give a gauge invariant result, but otherwise 
play no role in the calculation). Furthermore, we already 
know that = 0, and hence the value of G2 must 
be fixed in terms of Gi and G4 through Eq. (|3.16p . so 
we need not calculate it explicitly. This leaves only Gi 
and G4, whose values can be extracted from (|5.16p by 
separating out the terms dependent on g^'^ and q 



This leads to 



Gi = 

G4 = — —j= fv 

where the overlap integral is 



(5.17) 



Id3= bCk,q)^J^HP+.k)ll^fj{P-,k). (5.18) 
Jk 



From Eqs. ([STM]) and ((37T5)) we obtain 



f,lD3 (5.19) 



3V3 Af + m-' 
G^'iQ') = -^Tl^f^^Ds (5.20) 



3V3M + m-' 



(5.21) 



Although formally different from the integral involved 
in the S-state transition, it can be shown that the integral 
Xd3 goes with 1/Q* for large Q^. The proof follows the 
lines presented for case I of Appendix G from Ref . [11] . 
As consequence G^^ — G^^ ~ l/Q**- 



2. Nucleon(S) A(D1) 

For the Dl transition the p term is no longer zero, and 
using Eq. (|5.14p . the property of the S = 1/2 projection 
operator, and the definition of O2 gives 



P Jk { m 



(5.22) 



f2-I 



Dl 



which contributes only to G2 

2^/3M 
G2 = — 

p md- 

where the Dl overlap integral is 

b{k,q) ^fl(P+,fc)^^(P_,fc). 



(5.23) 



(5.24) 



Comparing this calculation with Eq. (|4.6p , and using the 
connection j2 — > G0/2-, gives an explicit expression for 



Pi 1 

2 2 



, 2M 



'-D1 



1 

30^ 



2Mm 



'-Dl- 



(5.25) 



Next, using Eq. (|5.15p for S = 1/2 case (where the 
r^a^p term does not vanish), the Oi term for the Dl 
transition is 



-V3/„ / {b{k,q) 



g2 



x\2g 



P^[2P^ -2M7''] 
mM + b 



{Vl/2)p'a 



(5.26) 



From this we must extract the contributions to Gi, G2, 
and G4 [again ignoring G3 which, using the modified cur- 
rent (|4.13p . will be given by the gauge invariant condi- 
tion]. It is easy to see that G4 = 0, and 



Gi = 



Dl 



2M 
\/3 d+ 
8Af(2m + M) 



Dl 



^/3d+d- 

These contributions combine with (|5.23p to give 

2 ? 



(5.27) 



fv 



Dl 



3V3 M + m 



GB\Q'' 



3^/3M- 
4toM 

JC-'-Dl 



Dl 



V3Q2 

where fc is the analogue of (|4.9p 



fc — fi- - 



2m{M 



(5.28) 



(5.29) 



Note that the expression (|5.28p for Gq^{Q^) is consis- 



tent with (|4.10p if we use the expression (|5.25p and the 
connection jc — > Cq fc- 

Finally, as we have already discussed, the possible sin- 
gularity in G^ must be canceled by imposing the require- 
ment 



lim 2]ji 



(5.30) 



where A is a constant. This constraint predicts that the 
Dl contributions to the magnetic and electric form fac- 
tors will be zero at ~ 0. However, the constant A 
in the limit (|5.30p will in general be nonzero, predicting 
that Gq is finite as — > 0. 

For large Q^, we can write Tdi as a difference of two in- 
tegrals of the type Tds with different coefficients. Hence 
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Idi goes like 1/Q^, which gives a l/Q"* behavior for 
G^^jG^^ and G^^ ^ (because fc constant as 

Q2 ^ oo). 

In the overall, the asymptotic expression for the form 
factors are consistent with pQCD 



C. Sum of all valence contributions 

Considering the sum of all valence quark contributions, 
we obtain the contribution of the quark core, which we 
denominate by 'bare' (B) contribution: 



G^ 



G 



D3 
M 



bG 



Dll 
M I 



NlaGf,' 



bG 



Dll 
M I 



Gg(Q2) ^ NbGS\ 



(5.31) 
(5.32) 
(5.33) 



where we used the relations between the electrical and 
magnetic components for each state. Note that there are 
only two contributions for Gf , and one of them (G^/) 
is zero for = 0. As for G^ there is only the Dl state 
contribution. 

For completeness we mention here that the nucleon 
could also have a D-state. However, the nucleon (with 
total angular momentum J = 1/2) can only have the D- 
state with core spin 3/2. This nucleon D-state can be 
build using the ideas presented in the previous sections 
and leads to an additional contribution to Gq. We have 
not considered such a D-state admixture in this paper 
because the nucleon form factors can be well described 
at low ^ without including it. 



Note that this scheme is model dependent, because 
the decomposition in background and resonances ampli- 
tudes is not unique 18, 19, 43, 99]. However, once estab- 
lished the pion production mechanism (ttNN amplitude) , 
we can split G* in two contributions in a given formalism. 

Although our main goal here is to study the D-state ef- 
fects in the core valence quark wave function, even a qual- 
itative estimate of the D-states effects requires a simula- 
tion of the pion cloud effects. An effective parametriza- 
tion of the pion cloud in G%j was already introduced in 
a previous work ^68] . For G^ and Gq we consider in 
the present work t he p arametrization introduced in ref- 
erences [2i,[5ll73,[71l, which we will sketch now. 



A. Pion cloud parametrization of Gc 

In a pure SU(6) model the neutron electric form factor 
Geu would be identically zero and the multipolcs E2 and 
C2 in the jN — > A transition negligible. In the real world 
Geu is small but non-zero. 

Considering a constituent quark model with a con- 
fining harmonic oscillator potential with also pion- and 
gluon- exchange between quarks, Buchmann [53,[l2| con- 
cluded that the Geu data can be explained considering 
a two-quark current, with a quark-antiquark pair inter- 
acting with the external photon. In this description the 
neutron spatial extension, expressed in term of its radius, 
can be written as 



(6.2) 



VI. PION CLOUD CONTRIBUTION TO THE 
FORM FACTORS 

The previous section presented the contribution for the 
form factors from the photon-quark interaction in rela- 
tivistic impulse approximation, and within the spectator 
theory. But the description of the electromagnetic A^A 
transition requires also the presence of non- valence de- 
grees of freedom, which may involve two-body currents 
and/or sea quark contributions - dominated by virtual 
pion states, the pion cloud effects. 

In the language of the dynamical models, where the 
hadronic interactions are described in terms of a baryon 
core which interacts with mesonic fields, a transition form 
factor can be separated into two terms [H, [H, [11] : the 
contribution of the quark core, or Bare contribution, and 
the contribution from the pion cloud: 



G*JQ')^G^iQ^) + GliQ'), 



(6.1) 



where a holds for M,E,C and G^ denotes the correspond- 
ing mechanisms involving at least one intermediate pion 
state. This contribution is related with the long range 
interaction, while G^ contains the short range physics 
[4^ parametrized by the baryon wave functions. The 
decomposition (|6.ip was also considered in Ref. (t^ . 



where bq is the quark core radius (oscillator parameter). 
For the experimental result: ~ —0.113 fm^, we can 
estimate sa 0.6 fm. Within the same formalism, one 
concludes ^3 that 



G^(0) = -J—Mm^. 
' V M 6 



As, for low Q^, we can write for Ge„ 



GsniQ^) - 
and we obtain, for small Q^: 



^2 ' n 



(6.3) 



(6.4) 



(6.5) 



The relation (|6.5p can alternatively also be constructed 
from relations between the nucleon and nucleon to A 
transition magnetic moment, in the large Nc limit [23 ]. 
for low Q2 (Q^ <:<:l GeV^). 

Following Buchmann again, from a different perspec- 
tive [t^ . the nucleon form factors can be described by 
a symmetric quark core distribution plus an asymmetric 
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FIG. 3: Comparing the GJ/ data with the parametrization of 

Eq. dnn). 



pion cloud around the inner core. Considering the proton 
electrical form factor in particular, we can write 

GEp{q') = Gep'{Q^) + Gl;p(Q2), (6.6) 

where Gep'{Q^) is the 'bare' proton charge form factor, 
and G'^p(Q^) is the contribution due to the pion cloud. 
In the same picture the neutron electric form factor is 
however just given by the pion cloud and we may write 



(6.7) 



since the charge distribution in the neutron 'bare' core 
is zero. In the = limit Eqs. (|6.6p and (|6.7p are 
directly related with the nucleon radii. From Eq. (|6.4p 
we obtain G£;„(0) « —0.113 fm^. As for Gep we 

may write Gep{Q^) — 1 — r^^, where is the proton 
electrical squared radius. Now, can be decomposed as 
— « 0.78 fm^, where r^, « 0.67 fm^ represents the 
radius of the bare proton, the size of the proton being 
increased by the pion cloud. 

As Geu is determined by pure pion cloud effects, we 



conclude that Gq (|6.5p is the result of pion cloud effects 
(or equivalently, the Coulomb quadrupole form factor in 
the A transition would be zero, for the case of no 

pion cloud effects). The previous derivation assumes no 
contribution from the inner core (symmetric distribution 
in the core). This assumption is not valid in general 
but can be a good approximation for a small D-state 
admixture. We will therefore use 



GUQ') 



12m Ge. 
-Mm 
M 



(6. 



to represent the contribution of the pion cloud for . 

To check the consistency of this assumption, and be- 
fore using it together with the bare model built here, we 
compare the Rsm data with the results extracted from 
the electrical form factor data using the parametriza- 
tion (|6.5p . To estimate G*j^_j{Q'^) at the respective mo- 
mentum Q^, we consider the simple phenomenological 
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FIG. 4: Comparing the Rsm data with the prediction of 
Eq. II6.5P using the neutron electrical form factor data of 
Ref. [j. The nucleon model line corresponds to the model 
of Ref. [I]- The nucleon wave function parameters are pre- 
sented in table m 



parametrization of Ref. [2l| 



G;,(Q2) ^3Gcexp(-0.21Q2)^ 



{M + m)2 



(6.9) 



The quality of this parametrization for G\,j is presented 
in Fig. [3l The results are presented in Fig. [H where 
we calculated Geu from our spectator constituent quark 
model. Although according to Eq. (|6.5p the pion cloud 
contribution to decreases with Q^, its effect is not 
observed in the figure, due to the kinematic factor 
present in Rsm- 

Due to the nature of the derivation of Eq. (|6.8p (large 
Nc limit and ~ 0) we cannot say for sure whether or 
not the discrepancies in Fig. [4] are the result of the crude 
estimation (0(1/7V^) correction to the large Nc limit) 
or the result of neglecting the bare quark contribution. 
Reference [t^ estimates the D-state effects from this one- 
body current to be 20% to the final result. 



B. Pion cloud parametrization of 

Considering the large Nc limit, Pascalutsa and Van- 
derhaghen [24| related G^ and G^ at the photon point 

{Q^ = 0) 



Gc(0) 



4A/2 



M2 



;GUO)- 



(6.10) 



Using the relation (16. Sp between G^ and Gbk, and ex- 
tending the results for finite Q^, one has [2^ : 



G*e{Q' 



m\3/2 M2 -m2 GEniQ^) 



M 



(6.11) 



This result was derived in Ref. [2J], in the = limit, 
and must be restricted to low (Q^ << 1 GeV^). The 
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comparison between the Rem data and the predictions 
from Eq. (|6.1ip using the Geu data is presented in Fig.O 
In the Gq case, the exact SU(6) symmetry would im- 
ply Geu = and there is no contribution for the elec- 
tric quadrupole. In contrast, one cannot conclude that 
Eq. (|6.1ip results uniquely from pure pion cloud effects. 
In the large Nc analysis, both G% and G^ are 0{1/N^), 
to be compared with = 0{N^), which is estimated 
in terms of the magnetic form factor of the neutron 
[65l Ft^ . In that limit the valence quark core is domi- 
nant, but the next order correction can be originated by 
pion cloud effects or by angular momentum excitation of 
a quark. But, because G*^ can be written in terms of rf^ 
(or GEniQ^) for ~ 0), we will take (|6.1ip as the pion 
cloud contribution for G*^ for low Q^, neglecting next 
order corrections 0{1/N^)) in the large Nc limit. 



3/2 M2-m2 G£„(Q2) 



2V2 



(6.12) 



Reference [SSi] estimates the contributions due to the 
quark-antiquark states to be 88% of G*^ for = 0. 

The relation between G*^ and G^ represented in 
Eq. (|6.10p is also known as the long wavelength limit 
for the ratio G^(0)/G^(0)). It is the result of the condi- 
tions: Q"^ « KP-m? « M^, rri^, like in large Nc limit, 
where Af — m ~ 0{1/Nc)- Equation (|6.10p is also used 
to relate the electrical and Coulomb bare quadrupoles in 
the SL HI and DMT (4^] model s. 

A direct consequence of (|6.10p . if G^ and G^ are the 
only contribution for the respective form factors, is that 

M 



Re Alio) = i?sM(0). 
Note that the pion cloud contributions 



(6.13) 



and (jO^ 

for Gq and G^ respectively, goes with l/Q^ for large Q^, 
competing with the bare contributions (1/Q^ and l/Q'^ 
respectively) . [Assuming as in Ref. thatGsn - 1/Q^]. 



I3\, P2 c+,c_ d+,d- X,ms/m 
0.049 4.16 -0.686 1.21 
0.717 1.16 -0.686 0.87 



TABLE L Parameters of the nucleon wave function (/3i,/32) 
and quark form factors corresponding to the model II of Ref. 
[3. In each case we kept k+ = 1.639 and k_ — 1.823 in order 
to reproduce the nucleon magnetic moments exactly. 



As a consequence, the pion cloud contribution does not 
change the asymptotic behavior derived for G J and Gg . 
We need to have in mind, however, that the results for 
G^ and GJ; are derived under the assumption that is 
small. Buchmann [7l|, It^I argues that nevertheless, the 
pion cloud description for G^ can be extended also to 
the intermediate region (Q^ ~ 4 GeV^). 

With the parametrization of the pion cloud mecha- 
nisms using the Equations (|6.8p and (|6.12p we preserve 
the covariance of our calculation because Geu is evalu- 
ated using a spectator model 



VII. RESULTS 

In this section we present the numerical results of our 
model to the 7iV — > A transition. For the quark current 
we adopted the quark form factors from Ref. [3| based on 
a vector dominance model (VDM) parametrization: 



/i±(Q') 

/2±(g') 



= A + 



(1-A) 

1 + Qymi 



(1 



r/Mif 



(l~d±) 



l + Q2/m2 1 + Q2/Af2 



(7.1) 



In these expressions and are the masses of the 
vectorial mesons. The lower mass, Wi, = nip (or m^^), 
describes the two pion resonance (three pion resonance) 
effect and Mh = 2m, takes account of all the larger mass 
resonances. The parameter A was adjusted to give the 
correct quark density number in deep inelastic scattering 
0, [6^ . All the other parameters are presented in Table 

m 

We will divide this section into two subsections. First 
we consider the effects of the valence quarks. In partic- 
ular we test whether the bare contributions alone calcu- 
lated as explained in Sec.|V]can describe the experimental 
data. In the second subsection we add the effects of the 
sea quarks (pion cloud effects), with the phenomenolog- 
ical, parameter free, description of the pion cloud pre- 
sented in the Sec. IVII 



A. Valence quark contributions only: Models 1-3 

With only valence quark degrees of freedom the iVA 
electromagnetic transition form factors are described by 
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and 2. G*m dat a from from CLAS/Jlab 
[, DESY hoi arrd SLAG [lol. Rem data from MAMI 
LEGS [ll, MIT-Bates [li] and Jlab [3, [H. The 
"bare" data for Gjv^, shown in the top panel, were extracted 
from the SL analysis, Ref. [13]. For the fit we doubled the 
bare data error bars shown in the figure to constrain Gfj, but 
the extra that results from the fit of Gfj to bare data is 
not included in any of the reported in this paper. 



Eqs. (|5.3ip - (|5.33p . The free parameters of our model are 
the admixture coefficients a, b and the momentum range 
parameters of the scalar wave functions (|2.40p - (|2.4ip . In 
a previous work we adjusted the S-state A wave function 
to the G\,j data considering also an effective pion cloud 
contribution as 



Gl,iQ^) = GZ{Q') + GlAQ'), 



(7.2) 



where Gf^ is the contribution of the quark core and G'^j 
the pion cloud effects, parametrized by 



{'SGd), 



(7.3) 



where = (l + (5^/0.71) is the nucleon dipole form 
factor, a cut-off and A^r a coefficient the defines the in- 
tensity of the pion cloud effect. The factor 3 was included 
for convenience: when = 0, G][,j(0)/G'^^(0) = A^r, 
then Att measures the fraction of pion cloud (G^^(O) ~ 3). 
The parametrization (|7.3p simulates the main features of 
the pion cloud mechanism: significant contribution for 



= Q- falloff with increasing Q^. For more details see 
Ref. [il. 

Here we extend the predictions to the subleading 
quadrupole form factors G*^ and expressed in the 
ratios Rem and Rsm defined respectively by Eqs. p.lOp - 
((3TT|) . We kept the parametrization ((7^ - (fr3| for G^^, 
however G£f is no longer determined only by the A S- 
state, but now also includes contributions of both of the 
D-states. For this reason the parameters originally fixed 
in the S-state fit are now readjusted. 

We considered the Gl, data from CLAS/Jlab [11,111, 
DESY [lOl and SLAG |loH. For the electromagnetic ra- 
tios Rem and Rsm we use the data from MAMI [l^, [IH , 
LEGS J^], MIT-Bates [H and Jlab [ll,[3. Although 
there is no inconsistency in the G^ data, there is some 
ambiguity in the Rem and Rsm data, dependent on the 
analysis. For the form factor information to be extracted 
one uses data for the pion photoproduction reaction cross 
sections. Those cross sections are interpreted in terms of 
an amplitude that includes a background and a resonant 
contribution. In the process, the extraction of the mul- 
tipoles depends on assumptions for the background and 
resonance parametrization. The multipolc resonant am- 
plitudes are then varied to fit the cross section data [3 ■ 
Kamalov et al. [H presented a re-analysis of the GLAS- 
2002 data [T3| with significant differences from the orig- 
inal data. Similarly the GLAS-2006 Rsm data [H for 
> 3 GeV^ shows a dependence on different from 
the recent MAID analysis [t^ of th e same data. The 
Rem analysis from Arndt et al. jl02j | is in contradiction 
with all the published results. More recently Stave (67| 
showed that there is a significant discrepancies in the 
extraction of E2 and G2 from the data using different 
reaction models like SL and DMT in the region < 1 
GeV^. [This discrepancy can be reduced by refitting the 
models within the range < 1 GeV^ only, which how- 
ever prevents the range of the application of the models 
for higher regions.] 

As a first step, model 1 fits only the G%f and Rem data 
(using, as in Ref. [6l|, the "bare" data extracted by the 
SL model [i^ to constrain the bare form factor Gf^). 
This fit (together with the fit from model 2 described 
below) is shown in Fig. [S] The Goulomb form factor 
predicted by model 1 (and not used in the fit) is shown 
in Fig. El 



The parameters that were adjusted during the fits are 
shown in Table |TT1 Although we did not fit the Rsm 
data, the coefficient b which determines the strength of 
the Dl state was adjusted during the fit. [As emphasized 
in the previous sections, only the Dl state can generate 
a non- vanishing Rsm, or Gq ^0.] As we see in Table HIl 
the best description of the G\j and Rem data requires 
a small admixture of the Dl state (0.2%). To check the 
sensitivity of the fit to the inclusion of the Dl state, we 
considered also another fit forcing b = 0. This defines 
model 2. As Table HIl shows, the admixture with 6 7^ 
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Model A^,A^ ai,Q2 Q3,Q4 Agi.Qs ^S,m xgm,Xrem Xrsm,X^ 



1 


0.450 


0.344 


0.1956 


1.025 


8.15% 


1.41 






1.46 


0.344 


0.1978 


0.1165 


0.17% 


4.39 


2.72 


2 


0.448 


0.350 






8.16% 


1.21 






1.53 


0.343 




0.0991 




4.90 


2.83 


3 


0.479 


0.343 


0.1567 


1.0087 


8.50% 


3.33 


11.84 




1.30 


0.350 


0.1574 


0.2218 


15.2% 


3.80 


5.45 


4 


0.441 


0.336 


0.1089 


1.0094 


0.88% 


1.41 


5.68 




1.53 


0.337 


0.1094 


0.1880 


4.36% 


0.99 


2.51 



TABLE II: Model 1 fits G*m and Rem- Model 2 fits the same quantities with b = (no Dl mixture). Model 3 fits all variables 
but restricts < 1.5 GeV^ for Rsm- Model 4 includes an effective pion cloud in both Rem and Rsm (for < 4.3 GeV^). 
All models also fit Gfj to the bare data (as shown in the figures) but the extra that results from the fit to bare data is not 
included in the reported in the last two columns. 
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FIG. 7: Rsm from models 1 and 3. Data from MAMI [10, [n 
LEGS 123, MIT-Bates [H and Jlab [li,[I|. 



improves the description of the data only shghtly (x^ of 
2.72 versus 2.83), meaning that the role of the Dl state 
is not decisive for G^^ and Rem- 

Figure [7] shows the prediction of model 1 for Rsm [the 
result for model 2 is zero]. We conclude that the Rsm 
prediction from model 1 is an order of magnitude smaller 
than the data. 

The next step is to try to fit the Rsm data as well, still 
using only the valence quark degrees of freedom. How- 
ever, because of the zero in fc, Eq. (|5.29p . we also predict 
a zero in G^, Eq. (|5.28p . Using the parameters of Ref. ,4], 
fc passes through zero around ~ 5.6 GeV^. This zero 
is completely at odds with the data. It is therefore im- 
possible to fit G"^ over the entire range, and at this 
stage we restrict the fit to Rsm to the low momentum 
region < 1.5 GeV^. This fit defines model 3. 

The results from model 3 are shown in Fig. [5] In the 
last panel of the figure we show Rsm for < 1.5 GeV^ 
only, the range used in the fit. Figure [7] compares the 
results for Rsm obtained from models 1 and 3 over the 
entire range. Note the unavoidable zero for model 3 
at ~ 5.6 GeV^. The first conclusion from model 3 is 
that the fit gives Rsm only within the region < 1.5 
GeV^, and that the fit is a poor one (high Xrsm)- Also, 



the quality of the description of the G^,^ data is affected, 
as we can conclude from Table |lTl by comparing Xgm 
obtained in model 3 with the corresponding values ob- 
tained in model 1 and 2. Note also that even the quali- 
tative description of Rsm provided by model 3 requires 
a abnormally large admixture of the Dl states (15.2%). 
All these observations show the intrinsic limitations of a 
pure constituent quark model. They can be overcome by 
adding pion cloud effects, as motivated in the discussion 
in the previous section. 

At low the data teU us that i?sA/(0) ~ -4%, which 
is equivalent to Gq(0) « 1.1 [this follows from p. lip 
with G%j{0) ~ 3]. Without a pion cloud such a result 
can be obtained in this formalism only by requiring a 
large Dl admixture. But even with a large Dl admixture, 
the valence quark contribution for Rsm changes sign at 
about ~ 5.6 GeV^, and we are led to conclude that 
the valence quark degrees of freedom are insufficient to 
explain the G^ data for large Q^- This conclusion is 
consistent with both constituent quark models and the 
results from dynamical models. 

Let us discuss now the numerical values of the range 
parameters a^. Remember that can be interpreted 
as a Yukawa range parameter [J, Q, with a smaller a 
parameter representing a larger spatial range. It is in- 
teresting to see that in all the models with D-state com- 
ponents, ai ~ a2- This finding differs from the results 
obtained previously for a A wave function with only an S- 
wave component [68|. For the pure S-state, ai « 0.3 and 
Q!2 ~ 0.4. Apparently, the introduction of the D-states 
sets a new long range scale, with a.^ ^ 0.1 to 0.2, for 
i = 3, 4, 5. These longer range scales are also nicely con- 
sistent with the notion that the D-waves are peripheral. 
Finally, it is worth mentioning that the similarity in the 
values for ai,a2 — 0.35, suggests that the S-state effects 
are somehow model independent. This feature suggests 
that, in the future, it might suffice to chose only 1 pa- 
rameter to describe the S-state of the A, showing that 
the improvement that accompanies the inclusion of the 
D-states is robust. 

Summarizing: a qualitative description of the G\.j and 
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FIG. 8: Model 3. Data from Figs. Hand Fig. [71 



can be obtained using a A wave function com- 
posed predominantly of an S-state with admixtures of D3 
and Dl states. A fit based on a quark core requires an 
8.2% admixture of the D3 state. The inclusion of the Dl 
state is not at all necessary to explain the G\,[ and G*^ 
data. To explain the data using only a quark core 
requires an unusual large admixture of the Dl state, and 
is only reasonable for low Q^, failing totally for > 2 
GeV^ (see Fig. [T]). The conclusion from models 1-3 
is that a quark core only is not sufficient to explain the 
~* A transition data, even when the A wave function 
includes admixtures of Dl and D3 states. 



B. A mixed description: valence quarks and a pion 
cloud: Model 4 

We now add a pion cloud contribution to the G*^ and 



form factors 



GUQ')^GS{Q') + GUQ'), 



(7.4) 



where the pion cloud contributions are taken from 
Eqs. (|6.12[) and (|6.8p . respectively. There are no ad- 
justable parameters in the pion cloud components. For 
the neutron electric form factor GEn{Q^) we use model 
II of the Ref. [i| (see table H}. The hmit of validity of 
the pion cloud formulas is restricted to low Q^, which 
led us to restrict our fit to Rsm to < 4.3 GeV^ re- 
gion. The bare contributions from valence quarks come 
from (|5.32p - (|5.33|) . Then, Gf is the result of the valence 
quark contribution involving the D3 and Dl states, and 
Gp comes only from the Dl state. 

The results for model 4, the fit obtained using the pion 
cloud, are shown in Fig. [9l with the parameters given in 
table (along with the results for all the other mod- 
els). The description of GJ, (and the corresponding ratio 
Rsm) by model 4 for > 2 GeV^ favors the recent 
MAID analysis [t^ over the original JLab analysis ^iB'l, 
and its success or failure will ultimately depend on which 
of these analysis survives further study. 

By comparing the x^s in table |TT] for models 3 and 4, 
one concludes that the inclusion of the pion cloud for G*^ 
and G*(j does indeed improve significantly the simultane- 
ous description of the data of these two more problem- 
atic observables. In model 4 the contribution of the pion 
cloud at Q2 = is 86.9% for G%, and 72.5% for G^. The 
most important D state is the Dl, with a small admix- 
ture of 4.4% compared to a tiny admixture of D3 state 
of only 0.9%. As for the dominant G^^ form factor, the 
pion cloud contribution estimated using an S-wave model 
does not change with the inclusion of D-waves (44.1% for 
model 4 to be compared with 46.4% from Ref. [6^). The 
magnitude of the cut-off, Aj, in the pion cloud contri- 
bution to G*j^ is decreased slightly (1.53 GeV^ versus 
1.22 GeV^). The addition of the pion clouds term to the 
D-wave states does not change the previously observed 
approximate equality ai ~ a2, although its reduces very 
slightly the value of these range parameters. 

Even though the pion cloud contributions dominate 
the description of the small form factors, our model sug- 
gests that the corrections coming from the valence quark 
sector are still important to obtain the best description 
of the data. This can be confirmed in Fig. [51 in par- 
ticular for Rsm- Leaving aside the discrepancy at high 
(which depends on the resolution of the differences 
between the recent MAID analysis and the older JLab 
analysis) the figure shows that the pion cloud contribu- 
tion, which in our work is parameter free, underestimates 
the data. 

While model 4 fits the overall data well, and is clearly 
the best model found, there are at least three ways in 
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FIG. 9: Model 4. Data from Figs. El Fig. [3 and the recent 
MAID analysis [75(1. 



which our theory and the experimental data could be 
improved. First, the validity of our model for the pion 
cloud can be questioned, particularly away from the pho- 
ton point, = 0. Second, the analysis of the data for 
< 0.15 GeV^ and for large is uncertain. If the low 
data is excluded from the fit, we found that a higher 
quality description of Rsm was possible, with ^ 1-2- 
That there is some legitimacy in excluding this data may 
be seen from a comparison between different data sets, 
shown in Fig. (TUl The experiments for = 0.121 
GeV^ from MIT-Bates and = 0.126, 0.127 GeV^ 
from MAMI suggest a large negative fraction for the 
Rsm (around -7 %), in contradiction with the data for 
= 0.2 GcV^ from MAMI, and also the very recent 
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FIG. 10; Rsm from model 4 with data and preliminary CLAS 
data using an Unitary Isobar Model (UIB), Sato and Lee 
model from [4l| (SL) and 43] (SL2). The preliminary CLAS 
data was not used in the fit. 



preliminary CLAS data Hlga,!! {Rsm - -5%). The 
final analysis of the CLAS data should clarify this point. 

Finally, our treatment of the valence quark sector can 
be questioned. The factor fc in Eq. (|5.29|) gives a zero 
in Gq, which implies a dominance of the pion cloud for 
^ 6 GeV^. To study the model dependence on this 
behavior of fc , we probed a change in the quark anoma- 
lous magnetic moments k±. We tried to suppress the 
large behavior of /2- by redefining k± — > k± Az+ga ; 
with A an adjustable cut-off. This reparameterization 
decreased the overall obtained for the description of 
the transition data, but also increased the for the fit to 
the nucleon form factor data. Clearly this effect deserves 
more study. 



VIII. COMPARISON WITH OTHER WORKS 

In general our results agree qualitatively with Refs. [H, 
pTOj and support the general idea that the quadrupole 
transition form factors are dominated by pion cloud ef- 
fectsjlpjpjil [el, [ei, [63| or quark-antiquark states 
[HIMlzllzIl- In particular our pion cloud contribution 
is consistent with both chiral perturbation and Lattice 
QCD estimations. Also, according to results of chiral 
perturbation theory Rsm logrriTr as which 

implies signifi cant pion cloud effects at the physical pion 
mass [13, iH, [io3|. Additionally, effective field theory 
calculations of Gail et al. [2l|, [tJ] predict a dominance 
of the pion cloud effects. The recent quenched and un- 
quenched lattice QCD calculations [6J] for pion masses 
TOtt > 0.35 GeV predict only a small fraction of the ex- 
perimental result for for low . This fact suggests 
as well that the pion cloud effects are dominant in G*q for 
low (the enhancement of the chiral loop corrections 
relative to lattice data for small pion masses was shown 
in Ref. 0). 

Constituent quark models, as the Isgur-Karl model 
[28l , include S-states and a D-state admixture of typ- 
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FIG. 11: Left side: Parametrization to the "bare" form factors from Sato and Lee (SL) [42|, Dubna-Mainz-Taipai (DMT) 
[SiEHl Dynamical Utrecht-Ohio (DUO) [46| ] models. The circles represents the experimental data from Figs. |6] and Fig. [7] 
Right side: Model 4 compared with SL and "bare" data. The ratios were evaluated using the parametrization of GJ,/ from 
(16.91). In both cases the "bare" data is from Ref. [431 . 



ically 1%, but predict only a fraction of the total G*^ 
and near = 0. This feature is also shared by sev- 
eral relativistic quark models [S^, [S^l and by the valence 
contribution of our model 4 (see Figs. [9] a nd pH) . The 
exception to this role is the work of Ref. [S^l where a 
manifestly Lorentz covariant chiral quark approach was 
considered. In that work the effect of the pion cloud is 
reduced to the order of 10%, which is compensated by a 
significant contribution of relativistic effects when com- 
pared with non-relativistic quark models. A discussion of 
the predictions of quark models for low can be found 
in Refs. [lillllei,!!^. 

We conclude this section by looking at the implica- 
tions of describing the pion cloud with a dynamical model 
(DM). Dynamical models assume that the complete elec- 
troproduction (or photoproduction) amplitude is the iter- 
ation of a kernel composed of the sum of bare resonance 
pole(s) in the s channel, plus "left-hand cuts" (arising 
from the angular average of t and u channel poles com- 
ing from the exchanges of mesons and baryons). Fit- 
ting the data with such a model fixes both the bare res- 
onance parameters and the parameters of the left-hand 
cuts which dynamically determine the background, and 



in this way allows one to extract bare "data" , or that 
part of the form factors that would be present even with- 
out the dressing produced by the rescattering of pions. 
Predictions of a pure CQM could be compared directly 
with this bare data, since both exclude the same physics 
- all pion rescattering mechanisms. 

The bare "data" extracted from the fits of Ref. [i^ are 
compared to various theoretical models in Fig. [TTJ Since 
the comparison involves electromagnetic ratios and not 
absolute quantities, we used (|7.2p to parametrize G*j^.j 
which is present in all the ratios shown in the figure. 
The left hand panels show the comparison of bare data 
to the parametrizations used by Sato and Lee (SL) [4l| . 
Dubna-Mainz-Taipai (DMT) [Ij, [H] and the Dynamical 
Utrecht-Ohio (DUO) ^ models. Note that there is a 
substantial difference between the "bare" data [i^ and 
the parametrization initially used in the SL model [4l| . 
particularly for Rem- 

The figure eloquently exhibits that "bare" contribu- 
tions are strongly model dependent, in their size and even 
their sign, differing from model to model. All that we can 
conclude from these results, considering also the observed 
experimental data, is that "bare" contributions and pion 
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= GeV^ 


G 1 G TP G"/"" / Cj 

tii 1 III o / o 


DMT [73] 


-5.7% -4.7% 


SL [73] 


33% 36% 


DUO* 


136% -42% 


Buchmann [53, 70] 


12% 20% 


Model 4 


13% 18% 


= 1 GeV^ 


G w 1 G T? GnlGn 


DMT* 


-8.8% -8.2% 


SL* 


56% 51% 


DUO* 


11% -15% 


Model 4 


17% 18% 



TABLE IIL Bare contribution in different models, estimated 
by subtraction of the pion cloud. In the lines labeled with an 
the total result was not available and model 4 was used. 



cloud contributions are both sizable (see Table HIH) . The 
exception is the DMT model, where the bare contribu- 
tions are almost negligible for low 5%). As for the 
SL model the bare contribution is 33% and 36% for 
E2 and C2 respectively at = [TJ] (corresponding to 
a pion cloud contribution of 77% for E2 and 74% for C2 
at the photon point). A compilation of the bare contri- 
bution for E2 and C2 for different models is presented in 
Table HIT] for (5^ = and 1 GeV^ . For a summary of the 
literature see also Refs. [gsI. ItsUtI] . 



IX. CONCLUSIONS 

In this work we introduce for the first time the D-states 
in the covariant spectator formalism for the description 
of baryons as a quark-diquark systems, and apply our for- 
malism to the description of the form factors of the elec- 
tromagnetic A'^A transition. Covariant formalisms pro- 
vide a correct treatment of boosts and rotations, which 
are important to describe correctly the kinematics and 
the dynamics in the intermediate region [Q^ ^ 4 
GeV^). There are two D-states for the A: one for the 
valence quark core of spin 3/2 (D3 state), the other for 
valence core of spin 1/2 (Dl state). We show that these 
D states have the correct spin structure in the baryon 
rest frame. Within this framework we show here that a 
consistent model, with orthogonal nucleon and A wave 
functions, predicts non- vanishing contributions for the 
Electric and Coulomb quadrupole form factors, an in- 
direct signature of the asymmetry of the valence quark 
distribution in space. 

However, we start by finding that the D-states con- 
tributions are not enough to explain the experimental 
data for G*^ and G^. An admixture of 8% D3-state 
can explain the Rem data, but the Rsu data cannot 
be explained without a Dl component. Importantly, al- 
though, is that even a very large admixture of the Dl- 
state cannot explain the high behavior of this ob- 



servable. This conclusion is consistent with results from 
other constituent quark models in the literature. 

With this established we had to turn our attention to 
the pion cloud effects. We find that the pion cloud con- 
tributions are essential to an accurate description of the 
7A'' A transition, and that our best model (model 4) 
gives a good overall description of the A tran- 

sition form factors. In this model we used pion cloud 
effects derived in the large Nc limit, containing no ad- 
justable parameters, and our fit predicts that the A wave 
function is the sum of a large S-state component with an 
admixture of 0.9% for the D3-state and reasonable 4.4% 
weight for the Dl-state. The pion cloud dominants G*^ 
and G^, with contributions of 87% and 73% respectively, 
at the photon point. Like the valence quark contribution, 
the pion cloud contribution is also covariant because it is 
based on a covariant description of the neutron electric 
form factor. As the pion cloud parametrization presented 
here can be justified only for low (Q^ < 1.5 GeV^, 
according with Ref. [2^). in the future we are planing 
to include an explicit relativistic calculation of the pion 
cloud, to replace the effective parametrization. 

The momentum distribution of the D3 state is deter- 
mined by one parameter (015 ~ 0.20), and that of the 
Dl state by two parameters which turn out to be nearly 
equal (aa ~ 0:4 ~ 0.10). These values are smaller than 
the one parameter (ai ~ q!2 ~ 0.35) required to repre- 
sent the S state, consistent with the picture that the D 
waves are more peripheral. 

We conclude with two notes which concern also future 
developments: 

1) We found that the quality of the description of the 
data is very sensitive in the regions < 0.2 GeV^ and 
at higher > 3 GeV^. Evidence for problems in the 
data come from an apparent inconsistency between the 
G*(j data for ^ o.l3 GeV^ and ~ 0.2 GeV^. The 
new CLAS data [H,!!!,!! for 0.16 GeV^ < < o.34 
GeV^ can be useful to clarify the situation. For high Q^, 
the understanding of the differences between the GLAS 
analysis and the MAID analysis will be also crucial for 
future progress. 

2) Presumably the most accurate estimate of the pion 
cloud effects comes from dynamical models which com- 
pute the dressing of the bare quark currents by pion 
rescattering to all orders. As a byproduct, these dynam- 
ical models determine the parameters of the valence, or 
undressed, quark contribution, which can be compared to 
a quark model without pion cloud effects. We observed 
that the results from the bare, pure valence quark form 
factors strongly depend on the pion production model 
(see SL, DMT or DUO models presented in the text). 
It becomes therefore crucial to use valence quark mod- 
els to estimate directly the bare form factors as func- 
tions of Q^, and to understand the nature of the valence 
quark distribution in the nucleon and A system, as we do 
here. Since the valence quark distribution dominates the 
largest 7iV ~* A transition form factor, G^ [11], dynam- 
ical models should use valence quark models as input, in- 
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stead of relying on phcnomcnological parametrizations. 
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APPENDIX A: SPIN STRUCTURE OF THE 
D-WAVE MATRIX ELEMENTS 

In this appendix we first show how the spin 3/2 func- 
tion V3/2 [defined in Eq. (|2.35p ] satisfies the special spin 
3/2 constraint 







(Al) 



and then prove the relations (|2.36p and explicitly con- 
struct the matrix elements in the expansions (|2.37p . 

To show that V3/2 satisfies the special spin 3/2 condi- 
tion, go to the rest frame P = (M, 0,0,0) and observe 
that 




A 

where the 2x2 operator is 

ax^x = (iAlA'||A,)fTv 

with 



"A, A \ 


XA 


; 






(A2) 



(A3) 



<7i£x = < 




(A4) 



A = -1 



Examining the As > cases shows that 
(as XA = crix+ = 

A 

Eh.) 



XA 



§'^0X4 



aiX-=0. (A5) 



Similar results hold for the As < cases. 

The orthogonality and normalization relations (|2.36p 
follow immediately from the orthogonality and normal- 
ization of the spinors and polarization vectors, and the 
unitarity of the CG coefficient matrix. The completeness 
relations follow from the normalization and orthogonal- 
ity relations, but it is instructive to prove them directly. 



To do this go to the rest frame and compute the matrix 
elements of the projectors in the spherical basis. Define 

^AV =£Ap(^s)a;3 4p (^6) 

and note that this matrix is hermitian and when mul- 
tiphed by the projection operator (M -I- f)/{2M) — 
i(l -f 7°) is of the block diagonal form 




(AT) 



where the 2x2 submatrices have the symmetry property 

CxpIx' (^x = P-x -y ■ (A8) 

Hence there are only three independent elements, which 
will be chosen to be poo, pn, and poi. Using the defini- 
tions of the projection operators (|2.25p . in the rest frame 
we obtain directly 
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1 

V2 




(A9) 



These same operators can be calculated from the def- 
initions (|2.35p and the expansions (|2.36p . In the rest 
system the spherical components of these expansions are 



J2el^Vso.iPAs)vUPAs)el-p 

Aa 

= ^(iAi lA|5As)(iA'i lA'I^As) 

As 

X ma(P,Ai)ua(P,A'i) (AlO) 
where the operator O has the form 



AA' 



"AA' 




(All) 



with the 2x2 matrix 

ofA' =E(5^i lA|5As)(iA'i lA'|5As)xA,xl; ■ (A12) 

As 

This operator has the same symmetry properties as pfy , 
and by explicit computation 
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(A13) 
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in agreement with Eq. (|A9[) . Similar agreement is ob- 
tained for the other matrices. 

We now turn to the computation of the matrix ele- 
ments in Eq. (j2.37p . Using the direct product represen- 
tations (|2.35p (and a similar one for Wj) 

Cd 2S = Vso. (P, a, {P, Aa ) 

= (i A 1 m|5A,) (i A 1 m'll Aa) A„rn' 

A 

= iV8^k2y„2^ ^ {2mi lm|lTO') 

A 

X (Im' iA||AA)(lm iAI^A,), (A14) 

where the CG coefficients guarantee that m = As — A, 
m' = Aa — A, and mg = m' — m = Aa — A^. The sum over 
three CG coefficients is eva luate d using Racah coefficients 
W . For the cases at hand If04ll: 



{21711, f TO|f to') (f to' i A|§ Aa) (f to i A|5As) 



= V3(25 + l)l^(2,f,|,i;f,^) (2to£5A,||Aa) 

,2'-z./, (Af5) 

V ^ 

giving Eq. (jOTl) . 



v2 



where z = cos 9 {9 is the angle between k and q) and 

k2 





afs _ Of/3 

f/3 

-2/3 



(B5) 
(B6) 



We can express (jB4p in a covariant form considering co- 
variant expressions for S3 and R"^ 

S3 b{k,q) 

i?"^ ^ R°'P{P^^P_), 

In particular we can write in the A rest frame 



b{k,q) 



3 (fc • q) 



~9 



R°''^{P+,P^) 



1^ 



2 

f 



(B7) 
(B8) 



The expression R"^ {P+, P-) is the only covariant ex- 
pression for R"^ compatible with (jB6[) , as can be showed 
considering the most general expression: 



APPENDIX B: INTEGRATION IN k 

When the currents associated to the A D-states are 
written there is a dependence in the tensor 

/"^(F+,P_)= / V'^f'{P+,k)^jrHP+,k)^UP~,k). 

Jk 

(Bf) 

The properties of I"^ in a Lorentz transformation follows 
the properties oi'D"'^: 



V°'P{P'^, k') = K''„KP v''P{P+,k). (B2) 



Then 



/"'^(P|, Pi) = A"^A'l/'"''(P+, P_) 



(B3) 



In these condition a covariant expression for I°'P{P+, P_) 
can be derived in a particular frame and the extended for 
an arbitrary frame using (jB3p . 

Consider Eq. (jBl[) in the A rest frame. In the A rest 
frame the scalar wave functions are independent of the 
variable ip. In these conditions we can perform the ana- 
lytical integration in ip, replacing the integral expression, 
by a equivalent integral with an integrand function inde- 
pendent of Lp. The result is 

^ j dpV"P{P+,k)^S3R'^^, (B4) 



R°'I\P+,P_ 



J pa pP 



ArP^P^ + A,P^P^ 

AP + + AP + - 



M2 



pa pp 



eg 



(B9) 



where a, h, c, d and e are functions of Q^. 

Similarly, the identity jBTj is the only possible covari- 
ant representation of S'3. Equivalent representation in- 
volving the factors q ■ k or q ■ k are reduced to q ■ k. By 
definition of k and q: k-q = q- k = q-k. 

As consequence of the representation (jB7|) and (|B8p . 
the integral (|B4p in the A rest frame can be represented 
by 



— / d^I?"'5(P+, k) - b{~k, q)R"^{P+,P^). (BIO) 



As (jBlOp is expressed in a covariant notation we can 
obtain the equivalent expression for a different coUinear 
frame considering an appropriate boost. 

Using Eq. ([BToI) we can write for any coUinear frame 



(p+ , p_ ) = p"'^ (p+ ,P-)Id (p+ , P- ) , 



— Da/3 



(Bll) 



where 



Id{P+,P-)^ / b{k,q)^^''{P+,k)i,^^{P^,k) 



(BI2) 
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